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1. Introduction 

We compute quark bilinear renormalisation constants (RCs), based on the ETMC Nf = 2 
dynamical quark action which consists of a tree-level improved Symanzik gauge action and twisted 
mass (tm) Wilson fermions at maximal twist (T)]. Our results are automatically improved in the 
spirit of ref. [^]. In section 2 we present a new (non-perturbative) method for the calculation of 
the scale independent RCs, Za and Zp/Z$, based on the use of two valence quark actions and a 
standard calculation of Zy within the tm valence quark sector. In section 3 we describe the RI- 
MOM calculation of all RCs (both scale dependent and scale independent ones). 

2. Calculation of the scale independent RC 

In this section we present a calculation of the scale independent RCs, namely Zy, Za, Zp/Z$. 
The evaluation of Zy is based on the PCAC Ward identity method (see refs. [^] for details). This 
calculation leads to very precise results. The computational method for Za and Zp/Z$ is new. It is 
based on the use of two regularisations for the valence quark actions. One is the standard twisted 
mass action, while the other is the Osterwalder-Seiler (OS) variant [Q]. In the so called physical 
basis these actions can be compactly written in the form: 

Sval =a 4 ^^(x)(7V-/7 5 rWcr + Hq)w{x) , (2-1) 

x 

with W cr = -§I M V* V M +M cr {r = 1), y = (u df , r = diag(r H r d ) and \l q = diag^ pL d ). The 
twisted mass case corresponds to r u = — = ±1, while the Osterwalder-Seiler case is obtained 
taking r u = = ±1. Sea quarks are regularized in the standard tm framework. 

Consider that, for the two different choises of the matrix r, we perform the following two 
axial transformations of the quark fields, namely (u,d) = exp[/(y5T37r/4)](w',</) and (u,d) = 
exp[i(y$7t/4)](u ,d), respectively. Each of the actions ( |2.1| ) transforms respectively into an ac- 
tion with the Wilson term in the standard form (no 75 and no T3). This is a rotation into the tm 
basis at maximal twist. However the tm action has a mass term of the form i/ixj/ , while 
the OS one has YsW r ■ Consider, now, an operator Or defined in the physical basis. Under the 
two axial trasformations this operator transforms into two operators, called Of and O^,, which, in 
general, are not of the same form. However the respective renormalised matrix elements between 
given physical states have to be equal up to 0(a 2 ) effects. This is due to the fact that in the con- 
tinuum limit each of them should coincide, up to 0(a 2 ), with the corresponding matrix element of 
the unique physical operator, Or- Therefore, if we call Zo t and Zq. the respective renormalisation 
constants for the two operators, we have: 

Z 0t {O f ) tm =Zo t {O t ) os + 0{a 2 ) . (2.2) 

Renormalisation constants are named, as usual, after the basis in which the Wilson term has its 
standard form. For maximal twist, the operator renormalization pattern in the physical and twisted 
bases is shown in Table 1 for both OS and tm formalisms. The primed operators refer to the tm 
basis while the unprimed ones to the physical basis and we have adopted the notation, Or = uYd, 
for both the primed and unprimed case. 
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OS case 




tm case 




{AR)n,ud 




— ZaA^ U(1 


{^R)ll,ud = ZvA-n^ud 


= -ttvV'^ud 






= Z V V n,ud 


(VR)n,ud = ZAV^ud 


= ~ iZ ^VMd 


{P R )ud 


= ZsP u d - 


= iz sS'ud 


{Pfdud = ZpPud z 


= ZpKd 



Table 1: Renormalization pattern of the bilinear quark operators for the OS and tm case at maximal twist. 

Calculation of Zp/Zs'. Our method is based on comparing the amplitude g K =< 0|P|7T >, 
computed both in tm and OS formalisms. We start by considering, in the physical basis, the correla- 
tor Cpp{t) = L x < uy=,d (x) dy=,u{0) >, which at large times behaves like Cpp{t) ~ [exp(— m n t) + 
exp(—m K (T —/))]. In the twisted basis, this corresponds to Cyy (t) in the OS case and C p i p i (t) in 
the tm one. Based on Table [T[ this translates into 

[ glt± ] C ° nt = Zp[ g ' K± r + 0{a 2 ) = Z s [g' K ] os + 0{a 2 ) , (2.3) 

from which the ratio Zp/Zs is extracted. 

Calculation of Za : We undertake the calculation of f % in both OS and tm regularisations. In 
the tm case we use the Ward identity evaluation of the decay constant: = 2pL q g 7l /m 1 l . Note 
that in this case no renormalisation constant is needed [Bt]. Thus the pion decay constant can be 
extracted from the large time asymptotic behaviour of Cpp{t) as it is discussed above. 

For the OS case we use the correlators Cpp and Ca p (with r u = r<i = ±1). The large time 
asymptotic behaviour of the former correlator has been discussed above, while the latter goes like 
Ca p(?) — 2^[ ex P( — m nt) — exp(— m n (T — f))]. Combining these, we can extract the bare OS 
estimate of the pion decay constant as f° s = %a p/ 'gn^in- Since the tm and OS determinations of 
the (properly normalized) decay constant satisfy the relation 

[h±] cont =f nn + 0(fl2) = Za/ OS + ^2) ^ 

an estimate of Za is readily obtained. Since all computations are performed at finite mass, the 
final results for Za and Zp/Zs are finally obtained by extrapolation to the chiral limit. Moreover, 
maximal twist ensures that cut-off effects are of order 0{a 2 ) ([§],[||]). 

2.1 Results 

Our configuration ensembles for Nf = 2 sea quarks have been generated at three values of the 
gauge coupling, j8 = 3.80, 3.90 and 4.05, corresponding to lattice spacings a ~ 0. 10,0.09 and 0.07 
fm. We have performed 240 measurements for the two smallest /3 -values and 150 measurements for 
the highest one. In order to significantly reduce autocorrelation times, correlators were computed 
every 20 trajectories (each having trajectory length equal to T = 1 /2). Five sea quark masses have 
been simulated at j3 = 3.90 and four at the other two couplings. The smallest sea quark mass 
corresponds to a pion of about 300 MeV and the higher one is just above half the strange quark 
mass. Eight valence quark masses were used at each coupling; the lowest ones are equal to the sea 
quark masses, whereas the others rise to the region of the strange quark mass. For the inversions in 
the valence sector we have made use of the stochastic method (one-end trick of ref. [|6|]) in order 
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J8 3.80 3.90 4.05 

Z A 0.72(2)(1) 0.76(1)(1) 0.76(1)(1) 

Z P /Z S 0.47(2)(1) 0.61(1)(1) 0.66(1)(1) 

Z y 0.5814(2)(2) 0.6104(2)(3) 0.645 1(2)(3) 

Table 2: The results for the scale independent RCs for three values of the gauge coupling. 



to increase the statistical information. Propagator sources are at randomly located timeslices. This 
turned out to be an optimal way to reduce the autocorrelation time. Typical plots on the quality of 
the signal for the RCs (for fixed values of the bare coupling and masses) are shown in Fig. [j] 
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(a) (b) 
Figure 1: Asymptotic behaviour of scale independent normalisation constants 



Three methods were implemented in the RC computation. The first consists in calculating the 
RCs at fixed value of the sea quark mass for a number of valence quark masses and taking the 
"valence chiral limit . Subsequently, the RCs were quadratically extrapolated to the sea quark 
chiral limit 2 . The second method consists in inverting the order of the two chiral limits. The third 
method is simply the extraction of the RCs from the subset of data satisfying \i va \ = fi sea , which 
allows to reach the chiral limit with one single extrapolation in the quark mass. Our results from 
all three methods are compatible within one standard deviation. We present preliminary data from 
the second method, which has fits of better quality, in Table ||. The first error is statistical while the 
second is systematic coming from the difference between the central values of the various methods. 
A final analysis will be presented in a forthcoming publication. 

3. RI-MOM calculation 

The RI-MOM method is a non-perturbative, mass independent, renormalisation scheme pro- 
posed in ref. [Q]. For a detailed presentation of various technical aspects see ref. [§]. In our case the 

1 First and second degree polynomial fits in \L va i have been performed. 

2 A quadratic dependence on a\l sea is expected from the form of the sea quark determinant, assuming that lattice 
artifacts on the RCs are not sensitive to spontaneous chiral symmetry breaking. However we have verified that a linear 
fit in n sea leads to compatible results. 
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scheme consists in fixing the Landau gauge and computing the momentum space Green function 

Gf( P ,p') = Y,{«(x)(uTd) d(y))e- i P* + V-y , (3 .i) 

x ,y 

for a general quark bilinear operator uTd (with r = A, V,S,P, T) and the propagator is written as, 

S q = Ex<?W?(0)>e _i/ "* with q = u,d. 

Then the forward amputated Green function, Af d = S u (p) Gf? (p, p)Sd(p)~ l , is projected by a 
suitable projector Pr (essentially a properly normalized Dirac matrix). The RCs, Zr and Z q are 
obtained by imposing the RI-MOM renormalization conditions 



Z^(Z tt Z,)- 1 /2 r ^(p)| p2=M2 =z¥(Z u Z d y^Tr\KfPr\\f =] ? = 1, Z q ±-Tr 



P 2 



1. 



q 12 

L ' J f-=\i 

(3.2) 

The computation is done for fixed quark masses. The results are extrapolated to the chiral limit. 
The renormalisation scale /I has to satisfy the condition: Aqcd <C H <C %ja. 

The RCs, calculated in the chiral limit in the way described above, are 0(a) improved at large 
momenta [§]. Moreover an analysis based on the symmetries of MtmLQCD and the 0(4) symmetry 
of the underlying continuum theory shows that Y^(p) and Y^(p) are separately 0(a) improved 
for all momenta. In order to increase the statistical information, we computed the following com- 
binations: Z r = (Zf + Zp) /2 and Z q = (Z u +Z d )/2. 

The scale dependent RCs (Zp,Z$ and Zj) are obtained at a reference scale \i§ = a~ l , by can- 
celling the scale dependence jJ., at a sufficiently high order in perturbation theory: 

Z r ( flj Uo) = ( Z T (apL)/C T (pL) ) Cr(jUo) • (3.3) 

Here Cr = exp f a (^ da [yf(a)/P(a)] and ft, j8 are the anomalous dimension of the operator and 
the beta function respectively. They are known at N 2 LO for Zj and N 3 LO for Zj and Zp [§]. 

It is known that the RI-MOM estimate of Zp is contaminated by the presence of a Goldstone 
pole In the twisted mass theory this problem also arises for Zs, though 0(a 2 ) suppressed. All 
these contaminations are removed in the subtracted Green function |(TT|l : 



rsub/ 2 „ „ \ _ MgJV(gJ%) ~ ^ 2 r P,s(p 2 ^q 2 ) sr.,, 

L P,S \P iV-qxiP-qi) — ~ — T \P A > 

F~qi H-q 2 

where jj, qi ,n qi are non-degenerate valence quark masses. 
3.1 Results 

The simulation parameters are the same as those of section [2.1[ The RCs, computed at fixed 
sea quark mass and several valence quark masses, are first linearly extrapolated to the valence 
chiral limit. Subsequently, the sea quark chiral limit is obtained by linear extrapolation in \l 2 ecv In 
Fig. U we show the effect of the Goldstone boson subtraction for Zp and Z$ for which the subtracted 



Green function of Eq. (3.4) has been used; we see that this has an important effect on Zp, while Z$ 
is almost unaffected, as expected. Moreover, we note from Fig. |2| that once the scale evolution has 
been perturbatively divided out, the scalar RC is indeed scale independent, while the pseudoscalar 
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/3 Za Zp/Zs Zy Zp Zs Zj 

3.80 0.75(3) 0.47(3) 0.62(4) 0.30(1) 0.64(2) 0.73(4) 

3.90 0.76(2)(1) 0.63(2)(3) 0.65(2)(3) 0.39(1)(2) 0.62(1)(5) 0.75(1)(2) 
4.05 0.77(1) 0.65(2) 0.67(1) 0.40(1) 0.61(1) 0.79(1) 



Table 3: RI-MOM results for the RCs. Z P , Z s and Z T are calculated at scale /J. Q = a' 1 (see Eq. (3.3)). The 
results at j3 = 3.80 and 4.05 are preliminary and the quoted errors are purely statistical in these cases. 

one is still subject to large discretization effects. These are removed by linear extrapolation, giving 
a Zp final estimate as the intercept of the fit. 

In Table |3| we show our preliminary results for the RCs; for j8 = 3.80 and 4.05 the results 
correspond to the lighter value of the sea quark mass only. For jS = 3.90 the results come from 
a full analysis in the valence and the sea sector. The first error is statistical and the second is 
systematic due to an estimate of the 0(a)-contribution to the quark propagator which induces an 
0{a 2 ) correction in the determination of the RCs. A better estimate of the systematic errors will be 
available once we finalize the analysis on all the three values of lattice spacing. 

A first comparison for j8 = 3.90 between the results of Tables || and ^ shows that the values of 
Za and Z p /Zs are in nice agreement and of comparable statistical accuracy. The corresponding Zy 
results, though compatible within the quoted errors, show that the PCAC Ward Identity estimate is 
statistically more precise 3 . 

We would like to note that combining the result of Z P /Z S from the first method with that 
of Zs from RI-MOM, an alternative evaluation of Zp can be obtained, in which the problem of 
the pseudoscalar Goldstone boson pole subtraction is avoided 4 . For example, for j8 = 3.90 this 
calculation gives Zp = 0.38 which is compatible, within the errors, with the corresponding value 

3 A slightly different determination based on the same WI taken between two one-pion states gives very similar 
results; for example, for ji = 3.90 it is found, Zy = 0.6109(2) jjl2J]. Moreover from the Table 1 we find that the value of 
the ratio (Z^/Zy ) 2 |^=4 05 is consistent with the one found in Jl 3[. 

4 We note in passing that Zp 1 = Zp is the quark mass renormalisation constant in the twisted mass theory. 
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given by the RI-MOM calculation (see Table 3). The results of a precise statistical analysis will be 

given in a forthcoming publication. 
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